Cosmological perturbations in Horava-Lifshitz theory without detailed balance 



Anzhong Wang 1 and Roy Maartens 2 
1 GCAP-CASPER, Physics Department, Baylor University, Waco, TX 76798-7316, USA 
2 Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth POl 3FX, UK 

(Dated: January 14, 2010) 

In the Horava-Lifshitz theory of quantum gravity, two conditions - detailed balance and pro- 
jectability - are usually assumed. The breaking of projectability simplifies the theory, but it leads 
to serious problems with the theory. The breaking of detailed balance leads to a more complicated 
form of the theory, but it appears to resolve some of the problems. Sotiriou, Visser and Weinfurtner 
formulated the most general theory of Horava-Lifshitz type without detailed balance. We compute 
the linear scalar perturbations of the FRW model in this form of HL theory. We show that the 
higher-order curvature terms in the action lead to a gravitational effective anisotropic stress on 
small scales. Specializing to a Minkowski background, we study the spin-0 scalar mode of the gravi- 
ton, using a gauge-invariant analysis, and find that it is stable in both the infrared and ultraviolet 
regimes for < £ < 2/3. However, in this parameter range the scalar mode is a ghost. 

PACS numbers: 04.60.-m; 98.80.Cq; 98.80.-k; 98.80.Bp 



I. INTRODUCTION 

Recently, Horava proposed a quantum gravity theory 
1], motivated by the Lifshitz theory in solid state physics 
2|. Horava-Lifshitz (HL) theory is non-relativistic 
and power-counting ultraviolet (UV)-renormalizable, and 
should recover general relativity in the infrared (IR) 
limit. The effective speed of light diverges in the UV 
regime, and this potentially resolves the horizon problem 
without invoking an inflationary scenario. In [3, 0| ; the 
general field equations were derived, and given explicitly 
for FRW cosmology, from which it can be seen that the 
spatial curvature is enhanced by higher-order curvature 
terms. This could open a new approach to the flatness 
problem and to a bouncing universe 0, @ . It was also 
shown that almost scale-invariant super-horizon curva- 
ture perturbations can be produced without inflation Q . 

Horava assumed two conditions - detailed balance and 
projectability. He also considered the case without de- 
tailed balance. So far most of the work [3|, U, la, |8H10|| 
on HL theory has abandoned projectability. However, 
breaking the projectability condition seems problematic 
[HKlIj]. With detailed balance, it was shown that mat- 
ter is not UV stable [5] . In addition, a non-zero negative 
cosmological constant is required, and it also breaks the 
parity in the purely gravitational sector . 

Because Lorentz invariance is broken in the UV, HL 
theory contains a reduced set of diffeomorphisms, and 
as a result, a spin-0 mode of the graviton appears. This 
mode is potentially dangerous and may cause strong cou- 
pling problems that prevent the recovery of general rela- 
tivity in the IR limit [l4Hl6|. 

In order to avoid these problems, one possibility is to 
keep the projectability condition. With this condition, 
Mukohyama argued that the problems found in pj], [l6| 
can be solved by the repulsive gravitational force due to 
the nonlinear higher curvature terms [ll| . In addition, 
without the pro ject ability condition, the theory seems to 
be inconsistent |17| . 



By abandoning detailed balance but still keeping the 
projectability condition, Sotiriou, Visser and Weinfurtner 
(SVW) showed that the most general such HL theory can 
be properly formulated with eight independent coupling 
constants, in addition to the Newton and cosmological 
ones [Hj]. Among these eight coupling constants, one is 
associated with the kinetic energy, which leads to the 
spin-0 scalar graviton, and the other seven are all related 
to the breaking of Lorentz invariance, which are highly 
suppressed by the Planck scale in the IR limit. 

In this paper, we study linear scalar perturbations of 
FRW models in the SVW set-up. The paper is organized 
as follows: In Sec. II, we give a brief introduction to the 
generalized HL theory formulated by SVW [l3j . In par- 
ticular, we add matter fields (not considered in [Hj]) and 
generalize the dynamical equations and the Hamiltonian 
and super-momentum constraints. Calcagni constructed 
the action for a scalar field @ , while Kiritsis and Kofinas 
considered the same problem, and then generalized it to 
a vector field [4|. However, the general coupling of mat- 
ter to this theory has not been worked out yet, since we 
no longer have the guide of Lorentz invariance. We shall 
not be concerned with this issue in the present paper, 
and simply assume that it can be done and represented 
by a general matter action, from which we can derive 
the conservation laws of the matter field. In Sec. Ill, 
we present the Friedmann-like field equations for FRW 
models with any curvature k. In Sec. IV, we first briefly 
discuss different gauge choices, and then study the lin- 
ear scalar perturbations of FRW models, working in the 
quasi-longitudinal gauge. We obtain the perturbations of 
the dynamical equations and the Hamiltonian and super- 
momentum constraints. We also compute the perturbed 
matter conservation equations. In Sec. V, we study the 
spin-0 scalar mode of the graviton in a Minkowski back- 
ground, by specializing the formulas developed in Sec. 
IV. We find that this scalar mode, which could poten- 
tially undermine the recovery of general relativity in the 
IR limit, is in fact stable in both the IR and UV regimes 
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for < £ < 2/3. However, it should be noted that in this 
range, the kinetic term of the scalar mode in the action 
has the wrong sign, as is evident in Horava's original re- 
sults pj, so that the mode is a ghost [HI [H, [3 ■ In Sec. 
VI, we restrict to perturbations of the flat FRW model, 
and find that the corresponding field equations are con- 
siderably simplified. In Sec. VII, we present conclusions. 



II. HORAVA-LIFSHITZ GRAVITY WITHOUT 
DETAILED BALANCE 

We give a very brief introduction to HL gravity without 
detailed balance, but with the projectability condition. 
(For further details, see [l3[ .) The dynamical variables 
are N, N l and g^ (i, j = 1, 2, 3), in terms of which the 
metric takes the ADM form, 

ds 2 = -N 2 dt 2 + g tJ (dx l + N l dt) (dxi + N j dt) . (2.1) 

The theory is invariant under the scalings 

t -> £ 3 t, x l -> lx l , 

N -> £- 2 N, N l -> r 2 N\ 9ij -> 9ij . (2.2) 

The projectability condition requires a homogeneous 
lapse function: 

N = N(t), N l = N*(t,x k ), g ij = g ij (t,x k ). (2.3) 

This is invariant under the gauge transformations, 

t = t + x°(t),x i = x i + X i (t,x k ). (2.4) 

The total action consists of kinetic, potential and mat- 
ter parts, 

S = ( 2 J dtd 3 xN^ (£ K -C v + C 2 C M ) , (2.5) 

where g — det <7y , and 

£ K = KijK* -(1-OK 2 , 

C v = 2A- R+^(g 2 R 2 + g 3 R ij R ij ) 



-j (//;/, ,:; + g 5 R /,',,/>"•' + g 6 R)R{R k 



CijR 1 ^ , where CV,- is the Cotton tensor. This term scales 
as £ 5 and explicitly violates parity. To restore parity, 
SVW excluded this term 

In the IR limit, all the high order curvature terms (with 
coefficients gi) drop out, and the total action reduces 
when £ = to the Einstcin-Hilbert action. 

Variation with respect to the lapse function N(t) yields 
the Hamiltonian constraint, 

J d 3 Xy/g{C K + C v ) = 8nG J d 3 x y /gj t , (2.8) 



where 



J = 2 [ N 



A I 



SN 



A I 



(2.9) 



Because of the projectability condition N = N(t), the 
Hamiltonian constraint takes a nonlocal integral form. If 
one relaxes projectability and allows N — N (t,x l ), then 
the corresponding variation with respect to N will yield 
a local super-Hamiltonian constraint Ck + Cv — SttGJ 1 . 
As argued in [l7j . this will result in an inconsistent the- 
ory. Results obtained by relaxing projectability should 
be treated with caution. 

Variation with respect to the shift N l yields the super- 
momentum constraint, 



8irGJ l 



(2.10) 



where the super-momentum tt' ij and matter current J 1 
are 



SCk 

Skj 



-K» + (!-£) Kg», 



,r = -N- 



.sc 



A I 



6N< 



(2.11) 



Varying with respect to gij, on the other hand, leads to 
the dynamical equations, 



1 



N^g 



(vVr y )' 



N k 



N K 7T 



+ -C K g v + F» 



-2(K 2 Y 3 

- 2-K k ^N^ 
+ 8ttGt 1 \ 



2 (1 - £) KK tj 



(2.12) 



+ - 4 [g 7 RV 2 R + g 8 (y iRjk ) (V*# fc )] . (2.6) where {R2 yi s R u K ^ s (/y + fjt) j% and 



Here C 2 = l/167rG, the covariant derivatives and Ricci 
and Riemann terms all refer to the three-metric g^ , and 
Kij is the extrinsic curvature, 



I Si 



V9 Sg 



'^J s =o 



(2.13) 



K. 



u - ^ i~9ij + Vi^ + Vj N t ) , (2.7) 



where iVj = gijN^ . The constants £, <7/ (I = 2, . . . 8) are 
coupling constants, and A is the cosmological constant. 
It should be noted that Horava included a cross term 



The constants are given by go = 2 AC, 2 , g\ = — 1, and 
n s = (2,0,-2,-2,-4,-4,-4,-4,-4). The stress 3- 
tensor is defined as 



2 S(y/gC M ) 
\[9 Sg i:j 



(2.14) 
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and the geometric 3-tensors (F s ) i - are denned as follows: that 



{F 2 ) l3 

ma 



- 2 9ij, 



R 



v 2 9w > 



2 {R ij - ViVj) i? - (R - 4V 2 ) i?, 
V 2 % - (ViVj - 3%) - 4 (i? 2 ) y . 

+ ^.9 iJ (3i? fc ;i? fe ' + V 2 i? - 2i? 2 ) , 

3 (iZy - ViVj) i? 2 - X - gij (R - 6V 2 ) i? 2 , 
(% + ViV,-) (#/w-R H ) + 2R (R 2 ) tj 

+V 2 (Jtffc,-) - V fc [Vi (Jlflj-k) + Vj (RRik)] 
-l 9ij [(R-2V 2 ) (i? WJ R«) 
-2V fc V, (RR kl )] , 



3 (i? 3 



3 r 



v 2 (n 



-v fc (v« (i? 2 ) jfc 
1 

"2^" 



RlRmRk' - 3VfcV; (i? 2 ) 



2ViVj (V 2 i?) - 2 (V 2 i?) % 
+ (ViiZ) (Vj-iZ) - ~<fe [(Vi?) 2 + 4V 4 i? 
V* (ViV 2 JZj + V,-V 2 iL fc ) 



V Rij - Vkyvi* T V j 

- (V<itf) (V,i?< ) - 2 (V fc i?0 (V fe %) 



1 



(VkRim) - 2 V fc V,V 2 i? fei . (2.15) 



The matter quantities (J*, J 1 , r IJ ) satisfy the conser- 
vation laws (Tij ]. 



<?«r fei - — (Vff J*) 
v9 



+ 



2Nk w 



N^g 



= 0, 



Ni 



v«r ik - — (v- 9Ji ) - -fv k r 
jk 

-— (VfciV, - ViN k ) = 0. 



(2.16) 



(2.17) 



III. COSMOLOGICAL BACKGROUND 



The homogeneous and isotropic universe is described 
by the FRW metric, ds 2 — -dt 2 + a 2 (t)jijdx I dx j where 

jij = (l + \kr 2 ^j 2 Sij, with k — 0,±1. For this metric, 
Kij = —a 2 H~fij and Rij = 2kjij, where H = a/a and 
an overbar denotes a background quantity. Then we find 



Ck = 3(3£-2)ff 2 , 

C v =2A- 6 4 + ^f + 2 -^f, (3.1) 
a z a 4 o° 

where f}\ = C^ 2 (3g 2 + 33) and (3 2 = 
C- 4 (9 54 + 3.g 5 +36)- 

Because of the spatial homogeneity, both Ck and £y 
are independent of the spatial coordinates, and the mat- 
ter quantities are 



J 



-2/9, J 1 = 0, Tij = pgtj, 



(3.2) 



where p and p are the total density and pressure. Then 
the Hamiltonian constraint (|2.8p reduces to the super- 
Hamiltonian constraint, Cjc(t) + Cvif) = 87rGJ'(t), 
which leads to the modified Friedmann equation, 



k 8ttG A 



— p+ 3 

2^/fc 2 4/? 2 fc 3 . . 

+ — + J V- ( 3 - 3 ) 



a* a 
From Eqs. (|2TTT|) and (|2TT3l we find that 



-A+ — 



fc 2/3ifc 2 12^ 2 fc E 



ff y = m-2)Hg ij . 



9°- 



(3.4) 



Then the dynamical equation (|2.12l) reduces to [ll 

, . a 8irG , 2 , 

(2-30- = — (p + 3p) + -A 

ad 3 

Afixk 2 16,3 2 fc 3 



(3.5) 



Similarly to general relativity, the super-momentum con- 
straint (|2.10[) is then satisfied identically, since J 1 = 



and, from Eq. (J33J), V^*- 



IT 13 I ■ 



0, where V» de- 



notes the covariant derivative with respect to 7^ . 

Using Eqs. (|3.3p and (13.51) . it follows that in the back- 
ground the matter satisfies the same conservation law as 
in general relativity, 



p + 3H (p + p) = 0. 



(3.6) 



This can be also obtained from Eq. (|2.16p . while 
Eq. (|2.17p is satisfied identically. 

In deriving Eq. (|3.3p we followed the usual assump- 
tion that the whole FRW universe is homogeneous and 
isotropic. In [llj |. it was argued that such an assumption 
might be too strong. If one relaxes the assumption and 
requires that only the observed patch of our universe is 
homogeneous and isotropic, one can introduce the no- 
tion of "dark matter as an integration constant" of the 
Hamiltonian constraint (|2.8p : pit) in Eqs. (|3. 31) and (|3 . 6[) 
can be replaced by pit) + £{t) in the observable patch, 
where £(t) = const/a 3 in the IR limit [HOii. Beyond 
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the observable patch, £ is necessarily inhomogeneous. In 
order to analyze perturbations on an FRW background, 
one needs to restrict the perturbations to the observable 
patch, which then raises issues about matching across the 
boundary of the observable patch. In our approach, the 
background is a homogeneous FRW spacetime, so that 
E = in the background. 



IV. COSMOLOGICAL PERTURBATIONS 

Linear perturbations of the metric give 

Sgtj = a 2 {q)h %j (f],x k ) , 

SN l = n i (r],x k ), 5N = 0(77)71(77), (4.1) 

where 77 is the conformal time. We decompose into scalar, 
vector and tensor modes 1201] . 



n = <j>, rii = B\i - Si, 
hij = -2ipjij + 2E\ij + Fi\j 



Fju+Hy, (4.2) 



Note that <f> is a function of 77 only, while B, Si, ip, E, F 
and are in general functions of both 77 and x k , with 
the constraints, 



s/ = o, f^ 1 = o, //; //, 

The perturbed energy quantity Eq. 

5 J 1 = -28fi. 



(4.3) 
is written as 
(4.4) 



In general relativity, Sfj, reduces to the density perturba- 
tion Sp. 

The perturbed matter current in Eq. (|2.11[) . on the 
other hand, decomposes as 



SJ* = \[q 



<1 



(4.5) 



and the perturbed stress tensor Eq. (|2.14j) decomposes as 
1 



or 



n 1 |( = o, n* = 0, ir j y = 



(SV + 2pip) 7 y ' + + 2U^ + U lj , 

(4.6) 



The angled brackets on indices define the trace-free part: 



(4.7) 



In general relativity, q 11 and q l reduce to the scalar and 
vector modes of the momentum perturbation —a(p + 
p)(v\ l + B\ l + v l - S l ), while SV reduces to the pressure 
perturbation Sp, and IT, II 1 and IP 5 reduce to the scalar, 
vector and tensor modes of the anisotropic pressure. 



A. Gauge Transformations 

Consider a gauge transformation as in Eq. (|2.4j) . with 

X°=e, X* + f| 4 = 0, (4.8) 

where C° = £°(r?), f = f (77, x k ) , £ = £ (r],x k ). Then 
the metric perturbations in Eq. (|4 . 2 p transform as 



B = B 



Si — Si + Fi — Fi — £j 



(4.9) 



where T~L = a' /a and a prime denotes d/drj. Note that 
these gauge transformations are precisely the standard 
forms given in GR. The only difference is that in the HL 
case, 4> an d £° are homogeneous. We can omit £° from 
B, since only the gradient of B occurs in the metric. 
However, we are free to maintain the £° term - and we 
do this in order that we can use the standard form of the 
gauge- invariant Bardeen potentials $, \& - see Eq. (|4.14[) 
below. Using the gauge freedom, we can restrict some of 
the quantities defined in Eq. 



Synchronous Gauge 

This gauge is defined by 

4> = 0, B = 0, Si = 0, 
and from Eqs. (|4.9[) we find that 



(4.10) 



f = l[ a <pd v +—,&= I Sidr) + d{x), 
a J a J 

t = J Bd 7] + J ^ (J aW) + C(x), (4.11) 

where C(x) and Ci(x) are arbitrary functions of x k with 

C i = 0, and Co is an arbitrary constant. Therefore, as 
in general relativity, this gauge does not completely fix 
all the gauge degrees of freedom. This gauge was used to 
study the scalar graviton mode in [l3j . 

Quasi-longitudinal Gauge 

In general relativitv, the longitudinal gauge is defined 
by B = E = F,i = [20]. However, due to the projectabil- 
ity condition, we see from Eq. (|4.9[) that we cannot set 
all 3 quantities to zero, although we are still free to set 
E = and Fj = 0. In addition, using the remaining de- 
gree of freedom, we can further set <fi = 0. Thus we can 
set 



0, E = 0, Fi = 0, 



(4.12) 
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with 

t; = - f acf>d V + ^, £ = E, £j = Fi, (4.13) 
fflj a 

which are unique up to a constant Co- We call this the 
quasi-longitudinal gauge (it has been used by [1, UU 
in the case where projectability is abandoned, and k — 0) . 

It should be noted that, as in general relativity, in 
each of these gauges only two scalars are left, and we can 
define the same gauge-invariant potentials as in general 
relativity [H 

$ = cP + -H{B-E') + (B-E')' , 

* = tp-H(B-E'). (4.14) 

(Note that in a different set of gauge-invariant vari- 
ables was used.) 



B. Scalar Perturbations 
in Quasi-longitudinal Gauge 

In the quasi-longitudinal gauge, the metric scalar per- 
turbations are given by 

ds 2 = a 2 [-dif + 2B\ t dx l dr] + (1 - 2tp) y ij dx i dx j ] . 

(4.15) 

Then from Eqs. flU) and (pTFJl . we find that 

K* = 
K = 



Kn + a [B\ij + (V/ + 2Hip) jij] 
K + a,- 1 (v 2 B + 3ip A 



= K lJ 



and 



B^ + -2Htp)j ij 
l K + 2HaT 2 (2 - 3£) (v 2 B + 3^' 

4ft & 



(4.16) 



C v = £v--t 1- 



+ 



3k) ^ 



48ftfc 2 



3k) ip 



24g 7 k 



C 4 a 6 



V 2 (V 2 + 3fcU. 



(4.17) 



To first-order the Hamiltonian constraint 

J d 3 Xy/j(5£ K + 5£ v ) = -16TTG J d 3 x^Sp. (4.18) 
Using Eqs. (|4T6| and (|4T7)) we find that 

(2 - 3£)U 



2k 



3k V 



2/3! 6/5 2 /c 3 57 ^2 



V 2 5 + 3i/>' 



(- 

V a 



C 4 a 4 



)(' 



3fc V 



4:TtGci 2 Sh =0. 



(4.19) 



The integrand is a generalization of the general relativity 
Poisson equation [2(| ■ Note that the Laplacian terms can 
be dropped from this equation, using the identity, 



d 3 x^V 2 f = 0. 



(4.20) 



At first-order the supermomentum constraint (12. 10)) is 



(2-3£)V/- 2kB-S,V 2 B 



87rGa qu 



(4.21) 



which generalizes the general relativity Oi constraint [20| . 
Note that in the general relativity limit (£ = 0) and in a 
Minkowski background (q = = k), Eq. (I4.2ip implies 



fl> = G(x), 



(4.22) 



where we used the Hamiltonian constraint (|4. 19|) to set a 
homogeneous function of integration to zero. This result 
is closely related to the fact that the spin-0 scalar mode 
of the graviton becomes stabilized in the limit £ = 0, as 
we show in the next section. 

The perturbed dynamical equations require the per- 
turbed (F s )y of Eq. (|2.15|) . The results are given by 
Eq. (jA.ll) in the Appendix. Using Eqs. (|A.2[) and (|4~6|) 
in Eq. (|2.12[) . we can find the perturbed dynamical equa- 
tions. The trace part gives 



tp" + 2Hip'-Fil) 
1 



3(2 - 30 ' 10 



V 2 B' + 2UV 2 B 



8nGa 2 
(2~ 30 



SV. (4.23) 



Here SFij 
Eq. (|A.1[) . and J- is defined as 



E9sC n 'S(Fsh, with S(F.) ij given by 



T = 



2a 2 



-A 



k 2ft/fc 2 12ftfc 3 



(2 - 30 
The trace-free part is 



B\ 



(ij) 



2UB\ 



(ij) 



SF, 



(ij) 



-87rGa 2 n 



(ij)- 



(4.24) 



(4.25) 



These two equations generalize the general relativity ij 
perturbed field equations (20j . 

The perturbed parts of the conservation laws (|2 . 1 6[) 
and (|2~T7l) give 



Vld 3 x 



V + 3U (SV + 5fi)-3 (p+p) $ 



q + 3Uq -aSV-^ (v 2 + 3fc) n 



= 0, 
(4.26) 
= 0. (4.27) 



The energy conservation equation is an integrated gen- 
eralization of the general relativity energy equation, and 
the momentum equation generalizes the general relativity 
momentum equation [20( | . 
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V. SCALAR GRAVITON ON MINKOWSKI 
BACKGROUND 

In the general relativity limit on a Minkowski back- 
ground, the scalar graviton mode should be suppressed. 
Otherwise the recovery of general relativity would be ob- 
structed. By contrast, when £ 7^ 0, we expect the scalar 
mode may play a significant role. 

We set a = 1, k = = A, J* = = J 1 and r ij = 0. 
Then 

SFij = - (l + aiV 2 + a 2 V 4 ) - %V 2 ^) , 

6C K = 0, SC V = -4V 2 t/., (5.1) 

where ay ee C~ 2 (8£ 2 + 3g 3 ), a 2 = -(~ 4 {8g7 - 3g 8 ). Note 
that ip = since H = 0, so that ■0 is gauge invariant. The 
Hamiltonian constraint (|4.18[) is satisfied identically. (It 
is interesting to note that if the projectability condition is 
given up, the Hamiltonian constraint becomes the super- 
Hamiltonian constraint SCk + 5Cv = 0, which gives the 
strong condition V 2 ip = 0.) 

The super-momentum constraint (|4.21[) gives 

(2 - 3£) ip - £V 2 B = 0. (5.2) 

Substituting Eq. (|5.1[) into the dynamical equations 
with / = 0, we find 

(2 - 3£) Utp + V 2 B) = 2 + ai V 2 + a 2 V 4 ) V 2 ^. 

(5.3) 

We consider the cases £ 7^ and £ = (general rela- 
tivity limit) separately. 

a. 

From Eqs. (15.21) and (I5.3[) we obtain the wave equation 
i> - 4 (l + "iV 2 + a 2 V 4 ) V 2 V = 0, (5.4) 
where c 2 ee £/(2 — 3£). In Fourier space, 

^„ + w 2 VVi = 0, w 2 ee n 2 c^ (l - ain 2 + a 2 n 4 ) , (5.5) 

where n is the wave-number. It is clear that the solution 
is stable in the IR limit, provided that < £ < 2/3, 
which is equivalent to 1/3 < A < 1, where A = 1 — £ 
is the parameter used in [l|. This is in contrast to the 
conclusions obtained in [2l| . in which it was found that 
ip is not stable for any choice of £. The main reason is 
that in [2l| the authors considered the case with detailed 
balance, or at most with 'soft' breaking (and not full 
breaking) of detailed balance. With the most general 
breaking of detailed balance [T3|, we find that to have 
stability in the IR, it is necessary that < £ < 2/3 (or 
1/3<A<1). 



In addition, ip is also stable in the UV regime for 
a 2 > 0. It is stable in intermediate regimes provided 
that either (a) a 2 > and a\ < or a\ < 4a 2 ; or 
(b) a 2 = and a\ < 0. 

B. £ = 

When £ = 0, from Eq. (|4~22|) we find that ip(t, x) = 
G(x). Inserting this into Eq. (|5.3[) and integrating, 

V 2 B = H{x) + 1 (l + ay V 2 + a 2 V 4 ) V 2 G{x), (5.6) 

where H(x) is an arbitrary integration function. Thus 
it appears that the scalar graviton has a growing mode 
oc t. However, it is important to note that B is not gauge- 
invariant and therefore B does not directly determine the 
stability of the spin-0 scalar graviton mode. The gauge- 
invariant variables defined by Eq. (|4.14p are in this case 

$ = B = I(x), * = ip = G{x). (5.7) 

Here I(x) is determined by Eq. (|5.6I) : in Fourier space, 
I n = (1 - otyn 2 + a 2 n A )G n . 

Clearly, neither of the gauge-invariant variables is 
growing with time. As a result, the spin-0 scalar graviton 
is indeed stable in the general relativity limit (£ = 0) on 
a Minkowski background. 

This conclusion appears to contradict the one obtained 
by SVW [13|. However, a closer analysis shows that in 
terms of gauge-invariant variables, the results are consis- 
tent. The synchronous gauge variables used in [l3j are ip 
and 

h = -6ip + V 2 £, B = 0, (5.8) 

and they find that 

E = L{x) + M(x)t + Q{x)t 2 . (5.9) 

The scalar mode appears to be growing because h is. 
However, by Eq. (|4.14l) . the gauge- invariant variables for 
the SVW solution are $ = -E = -2Q(x) and * = ip = 
G(x) - neither of which is growing. 

Our conclusion is also consistent with the results ob- 
tained recently by Mukohyama (Tlj . 

It is interesting to note that the coupling of the spin- 
scalar graviton to a dust fluid on a Minkowski back- 
ground does not alter this conclusion. In fact, one can 
show that ip and B will satisfy the same equations as 
above in both cases, £ 7^ and £ = 0. The only dif- 
ference is that now the Hamiltonian constraint (14. 19|) re- 
quires the matter energy quantity to satisfy the condition 
Jd 3 x5n = Q. 

VI. SCALAR PERTURBATIONS OF THE FLAT 
FRW MODEL 

We return now to an FRW background. In the flat 
case, k = 0, we find that the perturbation equations sim- 
plify considerably. 
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The super-momentum constraint (14.21| reduces to 



VII. CONCLUSIONS 



(2 - 3$)^' = CV 2 S + 8nGaq. 



.1 



Integrating it over space and using the Hamiltonian con- 
straint (I4.19j) , we find 



d 3 x (3Uq + aSfx) = 0. 



When k = 0, we also have 



} + + 3^ 4 ) - ^ 2 ) ^ (6 - 3) 

Then the trace-free dynamical equation (|4.25p gives 



(a 2 £?)' = (a 2 



aiV 



2 , a 2- 4 



V 4 V " 87rGa 4 n, (6.4) 



while the trace equation (|4.23[) reduces to 



V>" + 2Hip' 



2-3£ 



(1 + 5 



v 2 + ^v 4 ) vV 



87rGa 2 



35-p + (2 - 30v 2 n 



(6.5) 



3(2-30 

The conservation laws Eqs. (|4.26[) and (I4.26P reduce to 



d 3 5 



V + 3H (SV + Sfi)-3(p + p) ip' = 0, (6.6) 
2 (6.7) 



q + 'iUq = aSP + -aV z tt. 

o 



Note that not all the equations are independent. Equa- 
tion (16.51) can be derived from Eqs. (|6.1[) . (|6.4p and 
(|6.7p . Therefore, we are left with three first-order evo- 
lution equations, (I6.ip . (|6.4[) and (|6.7p . and two integral 
constraints, Eqs. (I6.2p and (I6.6p . for the six unknowns, 
tp, B, 5/j,, SV, q and IT. 

In terms of the gauge-invariant variables defined in 
Eq. (|4.14l) . we can rewrite Eq. (|6 .4[) as 



V s 4>. (6.8) 



The last term on the right acts as an effective anisotropic 
stress from HL gravity, i.e. from the higher-order curva- 
ture terms: 



IE 



1 



8nGa 4 



ai + ^iv 2 Wv 



(6.9) 



This stress is strongest on small scales, and is suppressed 
on large scales. It might provide a signal to distinguish 
the HL theory from general relativity. 



In this paper, we systematically studied the linear 
scalar perturbations of the FRW models in the SVW 
setup [ill , which is the most general theory of HL type [l[ 
when detailed balance is abandoned, but projectability is 
maintained (and so is parity). 

We generalized fl3j |. who considered only a vacuum 
Minkowski background. In addition to generalizing the 
geometrical terms, we included matter and derived the 
conservation laws. We have not specified the type of 
matter, as it is still an open question how to construct 
the matter Lagrangian Cm, although scalar and vector 
fields have recently been studied [4j, |5| . 

Working in the quasi-longitudinal gauge, we obtained 
explicitly the perturbed Hamiltonian constraint (|4.19l) . 
the super-momentum constraint (|4.2ip . and the dynam- 
ical equations (|4.23[) and (|4.25p . The perturbed conser- 
vation laws are given by Eqs. (|4.26p and (|4.26l) . 

A crucial issue in the HL theory and its generalizations 
is the spin-0 scalar graviton mode. By specializing the 
FRW background to its Minkowski limit, we showed via 
a gauge-invariant treatment that this mode is stable in 
the IR limit for < £ < 2/3. It is also stable in the UV 
regime, provided that the arbitrary coupling constants (77 
and gs are suitably chosen. The apparent contradiction 
with the results of [l3| is resolved via a gauge-invariant 
reformulation of their results. This is consistent with the 
results of Mukohyama [11]. We also showed that this 
conclusion is true when coupling the scalar graviton to 
a dust fluid in Minkowski spacetime. This result is also 
different from the one obtained in [2lj j . in which it was 
shown that the scalar mode is not stable for any given 
£. The main reason is that in [2l| the authors considered 
the case with detailed balance, or at most 'soft' breaking 
of detailed balance. 

The stability condition < £ < 2/3 has the unwanted 
consequence that the scalar mode is a ghost [il [llllill[l8l ]. 
To tackle this problem, one may consider the theory in 
the range £ < and then try to remove the instability of 
the scalar mode via the Vainshtein mechanism (22[. 

Our general formulas for the FRW background provide 
the basis for further work to analyze cosmological tensor 
perturbations, inflationary perturbations and large-scale 
structure formation in the framework of the generalized 
HL theory. We showed that there is an effective grav- 
itational contribution to the anisotropic stress on small 
scales, Eq. (|6.8p . so that in HL theory we have $ 7^ W 
even in the absence of matter anisotropic stresses. 
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Appendix A: Perturbed {F s ) i:j 

To first-order, the (i fl s ) i ■ are given by 



(Fo) 
(Fi) 



(F 5 



(Ft 



1 
2 

6k 2 



-a jij + a ^7y, 

3^|« + 7y (V 2 + 3fcj V 
^ (v.Vj - 7y V 2 ) (v 2 + 3*) V, 



2 7*7 

er a 



(Fsh = ^7, + ^(v 2 + 2fc)^ 

-^2 ( 3 ^ 2 + 10k ) (l>\ij ~ 7«VV) 



(%) - (vV) 



108fc 3 



a 

144fc 
36fc 3 
2k r 



7»i 



36fc 



- 3^|y + 7 y (5V 2 + 12fc) V 



V i V j -7i i V^) (V 2 +3fc)?/>, 
24fc 2 - 



+ 2jij (\7 2 + 3k) ^ 



3V 2 ViVj - 3 7ij V 2 ^ 



12/c 3 



-6 (^ViVj - 3 7lJ V 2 J V 2 ^ 

h2 7ii (v 2 - 37fc) VV 

12fc 2 . _ , , , . 

4fc ^7i 



6k 

'7a 



2(V 2 ^) r +V 2 



-7y V 2 (3V 2 +8fc) V 



|[(v 2 +3fc) 



3fc V 2 V> 



-7ij (v 2 + 3k) (v 2 + 2k) V 2 V>} 



(Fs) l3 = ^{^(V^O-V^V 2 ^) 

-V^-V 2 (V>| ife ) - 2ViVj V 2 (v 2 + 4fc) V 
+2 7 „ V 4 (v 2 +4fc) V 

+ 7lJ V fe V z V 2 (^,)}, (A.l) 
where V fc ?/> = i/ , ' fc ■ In addition, we also have 

+-j {(2 - 3£) 7 U + 2m/> - 2&/> 



B |ii -(l-C)7 ii V 2 J B }, 



(K 2 ) - (1 - f) = (3£ - 2) 



tf 2 



+2 



U(l-30.Bl y + (l-0 7 U 'V 2 B 
(2-30a(^-ff^>) f J '}, 
3(2 -30^ 2 M 



+ 



2{2 -ZCjH 



W Z B + 3a[^- Htjj 



7 U , 



7V fc V fe 7r 4i + 2n kii V k N : > ) + n ij V k N k 



(2 - 3£)H 



2Bl« - 7 ij V 2 B 



(A.2) 
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